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Abstract 


By examining the exact analytic solution of a kinetic model of collieion** 
al Interaction of ionospheric ions with atmospheric neutrals in the Bhatnagar-* 
Gross-Krook approximation, we show that the onset of intense auroral electric 
fields in the topside ionosphere can produce the following kinetic effects: 
(1) heat the bulk ionospheric ions to ~ 2 eV, thus driving them up to higher 
altitudes where they can be subjected to collisionless plasma processes; (2) 
produce a non-Maxwellian superthermal tail in the distribution function; and 
(3) cause the ion distribution function to be anisotropic with respect to the 
magnetic field with the perpendicular average thermal energy exceeding the 
parallel thermal energy. 
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I. INTRODUCTION 


Cround-based and satellite observations in the latter part of the past 
decade have brought about a major advance in the understanding of auroral arc 
formation processes: the electrodynamic interaction between the hot magneto- 

spheric plasma and the cold ionospheric plasma seems to play a central role 
[e*g«» reviews by Akasofu, 1981; Mozer, 1981; Kan and Lee, 1981; Chiu et al., 
1981]. Although a major consequence of this electrodynamic interaction is the 
production of a component of the electric field parallel to the magnetic field 
for the acceleration of electrons in discrete arc formation, a second, and 
possibly more far-reaching, consequence of the interaction is the finding ‘iliat 
the ionosphere is a significant source of plasmas in the magnetosphere as 
ionospheric ions are accelerated upwards by auroral electric fields related to 
substorms [e.g., Shelley et al., 1976; Mizera and Fennell, 1977; Richardson et 
al., 1981]. This is in addition to the polar wind [Banks and Holzer, 1968; 
Banks, 1979] which is a significant steady mechanism for transporting Iono- 
spheric plasma into cne magnetosphere. Actually, the idea that the ionosphere 
is an active participant in the magnetospheric response to solar-terreetrlal 
activity has had a fairly long history [e.g., Dungey, 1961; Axford and Hines, 
1961; Vasyliunas, 1970]. These authors pointed out the importance of the 
ionosphere in providing the appropriate Pedersen conductivity to close the 
magnetospheric convection circuit - a role wi.ich is crucial In modern theorlec 
of magnetosphere-ionosphere coupling [e.g., Chiu and Cornwall, 1980; Kan and 
Lee, 1980]. In connection with this role of the ionosphere in limiting mag- 
netospheric convection. Joule heating of the ionosphere [Walbridge, 1967; 
Fedder and Banks, 1972] and the neutral atmosphere [e.g., Ching and Chiu, 
1973; Straus and Schulz, 1976] have been considered not only as a thermal 
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energy source but eleo ee e eource to drive ionospheric motions which have an 
indirect Influence upon the magnetospheric convection flow. These studies do 
not address the question of the ionosphere as a subatorm*’related source of 
magnet oepheric plasma. In this paper we shall attempt to consider the kinetic 
properties of the ionospheric plasma as it responds to the onset of an 
enhanced auroral electric field, and as a source of magnetospheric plasma. 

The observational picture of the kinetic properties of such magneto- 
spheric plasmas of ionospheric origin is far from complete; thus, any serious 
theoretical effort at present must be in the category of "base building." 
Upward acceleration of ions in auroral electric and magnetic fields Is prch.ib- 
ly not difficult to understand since both the parallel electric field In 
inverted“V etructuree and the divergence of the magnetic flux tube favor 
adiabatic upward acceleration of ionospheric ions such as O'**. Nonndlabatlc 
features of auroral plaamas of ionospheric origin (such as heating and genera- 
tion of superthermal populations of ion beams and conics) are an entirely 
different matter. For ionospheric ions, the observations are particularly 
intriguing since not only are these ions somehow energised to superthermal 
energies in directions parallel [e.g., Richardson et al., 1981] and perpen- 
dicular [e.g., Klumpar, 1979] to the magnetic field, but the processes seem to 
operate over very wide ranges in energy (~ 6 eV - 10 keV) and in altitude (500 
km - 8000 km), and over a vide distribution of local times [Corney et al., 
1981]. Such preliminary observational results clearly indicate the direction 
of present and future theoretical studies of auroral ionospheric ions: how 
are superthermal and anisotropic ion populations formed and how are ions 
energised over some four to five ordert of magnitude in energy? 

Quite possibly, the answer to these questions may be in-aitu wave- 
particle interactions [e.g., Ungstrup et al., 1979; Okuda and Ashour-iU>dalla, 
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1981]; perhaps future quantitative alaulatione with nore realistic conditione 
will answer the question of foraation of beans (parallel energisation) and 
conics (perpendicular energisation) in the observed energies covering the 
ranges of 6 eV to > 10 keV. Whether such theories are realistic or not is 
outside the concern and scope of this par^r* Ra<^her» we are more interested 
in the orig in of such super thermal populations at the lowest energy ranges (< 
6 eV). The question of how special nonadiabatic features are formed out of 
the cold ionospheric ion population is quite pussling when one examines the 
conditions of the ionosphere. First, high-latitude ionospheric temperatures 
at F-region heights, as measured by the S3-3 satellite, are generally < 2500 '’ 
K ('>'0.2 eV) (Rich et al., 1979), far less in energy than the supertherraal 
fluxes observed. Second, Klndel and Kennel [1971] concluded that the unstable 
regions of electrostatic ion cyclotron waves were above the F-maximum and 
generally intlM far topside ionosphere ab«?ve 1000 km; this makes the relative 
abundance of 0*^ in such nonadiabatic populations (conics and beams) 
[Ghielmettl et al., 1978] even more pussling since 0*^ must be driven up to 
such altitudes by some preheating process so that 0'^ can at times be the 
dominant ion at altitudes > 1000 km* Note that in this situation, the abun- 
dance of 0^ at altitudes > 1000 km cannot be attributed to escape of H”*" to 
higher altitudes, as O'*" is also observed to be dominant at the equatorial 
regions during storm time [balslgdr et al., 1980]. Could it be possible for 
some kind of pre-heating process to operate in the ionosphere to drive up an 
intense super thermal (~ several eV) population of ionospheric ions into 
regions where ion cyclotron waves can act to energise these ions up to tens of 
keV? Recent observations (Locksood and Tltheridge, 1981] support such 
hypotheses. 
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As wc hsvs diacusssd hlchsrto, the probable occurrence of Ion Joule 
heating in the ionosphere by convection electric fields mapped doun to iono*- 
spheric heights [Fedder and Banks » 1972] has been well accepted. Since such 
calculations use a fluid approach, which automatically assumes a Maxwellian 
form for the distribution function, they cannot tell us about kinetic features 
of the ion distribution such as pitch angle anisotropy and superthermal popu- 
lations without going into extremely complex calculations with higher moments 
[Schunk, 1975]. In view of the necessity to understand how magnetospherlc 
ions may originate from the ionosphere and in view of the predominance of the 
ataady convection electric field as an ionospheric heat source, we are 
prompted to ask if the auroral electric field may not be the source of pre- 
heating «rhich provides the topside with superthermal ions. Since the observed 
nonadiabatlc features are kinetic in character, we are then driven to consider 
the kinetic reaponae of ionospheric ions to the onset of steady auroral elec- 
tric fields in a simple collisional kinetic plasma model - the Krook model. 
Analytically solubla models always overeimpllf y ; thus our main purpose is to 
gain soma insight into kinetic properties of the ion distribution under 
auroral ionospheric conditions rather than to attempt to "explain** the 
observed nonadiabatlc properties alluded to above. Our tiork, simple though It 
may be, differs from previous work on ionospheric heating using a fluid 
approach [Fedder and Banks, 1972] in that we are required to examine the non- 
Maxwellian features of the distribution function rather than to assume 
Hax%#ellian distributions at all times, as is done in fluid models. 
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II. FORMULATION AND PRELIMINARY CONSIDERATIONS 


Consider a uniform auroral electric field Imposed upon a horlcontally 
uniform Ionosphere consisting of Ions, electrons and neutrals whose distribu- 
tion functions are respectively: f^, f^ and f^. In the auroral Ionospheric 

region of Interest, where we assume Che magnetic field to be vertical and 
uniform (fi Bs), these three species Interact through collisions and the 
total self-consistent electric field The Interactions of these species are 
symbolically formulated in terms of three coupled Boltzmann equations for the 
distribution function f^^* k"(l,e,n): 

A 

where Ly Is the Boltzmann operator for species k of charge q|,e and mass mj^, 

+ V • 7 + (Fj^ a^) • 7^ E D/Dt (2) 

The force Fj^ consists of a combination of electric, magnetic and gravitational 
components. 


+ + e qj^(v * B z)/c 


(3) 


K. , is the appropriate binary collision operator. The self-consistent elec- 

^ .9 

eric field must satisfy Poisson's equation. 


7 • 8 - 4ve / d^v (fj - f^) 


( 4 ) 
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It It w*ll-known that (1) - (4) are very difficult to solve even vith the 

A 

sinplest assumptions for K. In the next section we shall show that a trun** 
cated form of these equations with a Krook-type collision cern can be solved 
to «ive some insight into the kinetic properties of ionospheric ions under the 
influence of a strong 

The highly symbolic formulation (1) - (4) nonetheless allows us to 
discuss certain asymptotic constraints relevant to the problem. If we iguore 
the fluctuation part of f (the auroral DC is usually observed to be much 
larger than the AC electric field), the drivers of non-adiabatlc effects are 

A 

the non-linear collision operators which allow energy to be exchanged 
among the three species. Thus, given sufficient time and no loss of particle 
species from the system, it is believed that the collisloual evolution of the 
system tends asymptotically toward an isotropic Maxwellian form for the 
distributions of all the species, although the equilibrium thermal energy of 
each species need not be the same. This means that the fluid moments 
approach, with the underlying assumption of Isotropic Maxwellian distribu- 
tions, is valid at times long after onset of E^; for ions the time scal.^; is 
likely to be sat by the ion-neutral collision time. In the E-region, where 
ionic Joule heating effects have been thoroughly studied [Fedder and Banks, 
1972] using the fluid equationc, the equivalent time scale as determined by 
the ion-neutral collision time is much less than one second but the neutral 
response time is 1/2 hour; therefore, we do not expect kinetic features such 
as superthermy and pitch-angle anisotropy to persist for much longer than the 
ion-neutral collision time after onset. The situation for the F-region and 
the topside (> 35C km) is quite differed.:. In this regime, the ion-neutral 
collision time (> 10 sec) is much longer than the ion gyration period (" 1/30 
sec); therefore, the expected evolution to Maxwellian isotropy takes place 


6 


over oany gyration perloda» perreiccing kinetic features to persist for tens of 
seconds or minutes* Can a strong drive superthermy and pitch-angle 
anisotropy In the Ionospheric Ion population during the evolution period 
between onset and the asymptotic Maxwellian state? It la the Intent of 
this paper to demonstrate that this question can be answered In the affirm- 
ative for the Krook model of weakly ionized plasma kinetics. 


Before we proceed with the solution of the Ionic segment of a Krook-model 
plaama in uniform electric and magnetic fields, it may be convenient to give 
an elementary discussion of the origin of the expected kinetic features 
supertheray and pitch-angle anisotropy. An Imposed causes ions to drift 
and collide with neutral at<»s of the thermosphere which act to deflect the 
uniform drift motion into somewhat random motion. In the initial stages (far 
from the asymptotic state) the rat** of change of Ion kinetic energy is not 
isotropic because the mobilities (iC|, of Ions In the F-reglon are 
different for action parallel and perpendicular to the magnetic field. The 
perpendicular energy changes according to [Alfven and Falthammar, 1963] 


- T, 


(5) 


where m. x, • v. /(v^ + Q^) (Rlshbeth and Garrlott, 1969] and Yi = 2/3 for the 
1 X In in 1 

F-reglon. In (5) and refer to mass and energy of the neutrals; Is 
the Ion-neutral collision frequency and Q la the Ion gyrof requency. The first 
term on the right hand sice of (3) Is due to the driving of the elec*’"^*: field 
and the second term la due to energy loss to the neutrals. The same amount as 
the energy loss In (5) appears In the corresponding equation for as neutral 
energy gain* The simultaneous solution of (3) together with the equivalent 
equation for does not concern us here. The point that the Ion population 
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may be anlsotroplc soon after onset Is made by noting that the parallel Ion 
energy is unlikely to Increase until sufficient perpendicular drift energy has 
been scattered through large angles into the parallel direction - probably 
after many collision times; whereas, by virtue of (5), increases in one 
or a few collision times. 

Obviously, the ideas discussed in this section are qualitative. For the 
rest of this paper, we shall consider the results of a simple solution of an 
ion component of the Krook model for weakly ionized plasma in steady electric 
and magnetic fields. 


111. ^tOOEL OF KINETIC RESPONSE 


The full tr««Cm«nc of th« kintcic retponto of « gravltatlonMlly" 
•tratlfUd wuakly-lonliod «od lnho»os«n«ou*ly-««gn*tlMd platat (th« F-r«glon 
lonosph«r«) to tha onaat of • strong electric field, ss fornslly presented In 
the previous Section, la exceedingly difficult to solve. It ts the purpose of 
this Section to Initiate a kinetic analysis of this t» ?nsitlon region between 
the colllslon'^onlnated E-roglon and the collisionless magnetosphere with <m\ 
oversimplified but exactly soluble version of (1) - (4). The slmpllf Icatii'us 
Introduced are] 

1. The neutral exponent !• assumed to remain In static equilibrium. 
Thus, the neutrals act as acatterers of Ion motion but do not pick up any 
energy In the process. This is a reasonable approximation if the Ionic con- 
centration Is sufficiently low or If the average neutral mass Is high. 
Neither condition Is strictly valid for the F-reglon. Fodder and Banks [1972] 
showed that in a fluid model the motion of the neutrals Is »n important 
determinant of the Ion heating on time scale greater than about one half hour 
after electric field onset; since we are dealing with ic.i response at early 
times (several collision times after onset), the assumption of iimnoblle 
neutrals Is approximately valid. We hope to relax this restriction In our 
next stage of kinetic model development. 

2. We Ignore the colllslonal Influence of electrons upon the dynamic 
response of Ions to steady electric field onset. This Is not a bad 
approximation for early tlsMS because Ion-electron collisions do not change 
the Ion energy by very much. 
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3. W« Ignore the fluctuating (AC) electric field in (A), The AC field 
driven by auroral electrona may be very important in ion heating, especially 
in the form of resonance with electrostatic ion cyclotron waves. This effect 
has been pointed out by others (e.g., Okuda and Ashour~Abdalla, 19tU). Here 
we are not interested in cyclotron heating but, as we have stated, concentrate 
primarily on how the bulk of auroral ionospheric ions at ~ 0.2 eV may be 
heated to several eV on the topside. A complete theory of ion conics and 
beams cannot, hotrever, ignore the effects of the AC electric field. 

A* The colllaion operator for the ion component of (1) is approKinwtU'd 
by the simple Krook model [Bhatnagar et al., 1934], which was ptoposed speci- 
fically to study the approach to equilibrium of weakly-ionised cullistonal 

plasmas. 


^in ■ V (x) (- f (x,v,t) + n (x,t) f^^ (x,v,t)l (b) 

where v is the ion-neutral collision frequency which thermalises the ion 
distribution function f towards the assumed isotropic Maxwellian form 
(Bhatnagar et a> , 19S4j, 

fy 5 (m/2» T(x,t))^^^ exp - (m/2T(x,t)](v - u(x,t))' 1?) 

Note that / d’*v f^ - 1. In (b) and (7) n, u and T are the density, flow and 

energy moments: 

n 8 f f (g) 

u s J d^v V f/n 


lu 


(9) 


( 10 ) 


T = / d^v mlv “ u)^ f/3n 

In order for the Krook model to coneerve the above three quantities, (6) - 
(10) have to be self'-consiatently Included in the solution of (1) [Bhatnagar 
et al., 1954J. The ion-neutral collision frequency v can depend on n, but we 
shall assume v to be a given function of z only* By using the Krook model, it 
is assumed that v does not depend on the velocity v because the krook model 
does not conserve number and energy if v • v (v, x). 

S. The plasoM is assuised uniform in the horizontal direction and g * 0 
in (3). Obviously, this simplification ignores auroral spatial scales, but 
little progress can be made otherwise* The distribution function f can thus 
depend only on z, the vertical coordinate. 

6* The driving DC field ^ is assumed uniform for t > 0 

t - tj^(t) - Ej^ e(t) X Ul) 

Since 7 • £ ■ 0, (4) implies / d v f ■ Oq where oq is the ionospheric electron 
density - assuised constant* Thus, self-consistency of the electric field (4) 
and assumption (11) implies n(x,t) ■ uq* This constraint can be v«>rlfied from 
our solutions* 

Our model of kinetic response is defined by (1) - (4) and (6) - (11)* 
With the minor exception of restrictions on S and v, thl? BK>del is essentially 
the extension of the Krook model to the case of collisional plasma in uniform 
electric and magnetic fields* We seek solutions of this model as an initial 
value problem in which Ij^(t) causes f to evolve with time from an initial 
isotropic Maxwellian state 
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( 12 ) 


f(x, V, t ■ Oi - Oq lm/2x exp - Iib/ZTq] 

where Tq is the cold Ionospheric temperature ("* 2500* K): T (x,0) ■ Tq, 

J(x. 0) - 0. 

Since the assumed force f is divergenceless in space, we seek a solution 
of our model by a phase'-space transformation (x, v, t) ♦ (x', V, t) so that 
D/Dt in (2) is transformed into 3/3t. This is acccjplished by the following 
transformation between the components, labeled (1, 2, 3), of the above 

vectors: 





'cos at - sin at 


'^ 1 ' 


- sin at 

1 

- 



+ (cE^/B) 


V* 

m ^ 


sin at cos at 

m 


.^2. 


cos at - 1 


X* - X3 - V3 t 


(13) 


(lA) 


(15) 




vl(T) 

vJ(T) 


e(T) 


(lb) 


Note that the coordinates (x* , v' ) are the time-reversed evolution of the 
coordinates of an ion under and Be starting at (x, v) at t * 0; thus 


rv'(x, V, 0)1 . 


Lx'(*. V, 0)J 

uJ 


(17) 


Further, since the transformation (14) between velocity spaces entails a 
simple rotation, the volume element is invariant: d^v ■ d^v*. The transla- 

tion in (14) amounts to a shift of origin in velocity space. Under this 
transformation, the Krook equation (from (1), (2) and (6)] 
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( 18 ) 


Of 

W 


0k 



become! (see Appendix A) 


3f(x',v',t)/at - v[*(x',v',t)l{nQfQlx(x'»v’,t), v(x’ ,v' - f(x',v*,t)}. (19) 


In (19), for the eeke of explicitness, we have written x, v as functions of 
(x’, v', t). Since we shell be working with the (x*, v’, t) coordinates we 
shall hereafter use x and v to denote functions of (x*, v', t) specified by 
the inverse transformation of (13) ~ (lb). 

The solution to (19), with the initial conditions (12), is 


f(J',v',t) - f(x',^,0) e“''*^ 


“o® 


-vt 


dTve 


vt 


fQ[x(x* ,v' ,t), v(x',v't),t1 (20) 


where v i v[s(x* ,v* ,t)] in general. Note that this integral representation 
solves (19) generally once f(x',v',0) and fQ(x,v,t) are given. General inte- 
gral representations for u(x',t) and T(x',t) are obtained by substitution of 
(20) into (9) and (10). Bacause of assumption 6 (l.e., 7*^^^ - 0 and (A)), the 
possible forms chosen for fo*^c constrained by 


n(x',t) - / d^v* f(x',v',t) - Uq 


( 21 ) 


We show in Appendix B that our choices (7) and (12) satisfy this constraint 
(at least for the simple model given below). The integral representation (20) 
shows explicitly that, as the consequence of collision with neutrals, the 
initial distribution f(x',v',0) disappears in a collision tiiM (1/v) while a 
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n«w convoluted dletributlon take* its place 


We have, in principle, solved the problen with (20) end (21); however, 
devising e tractable model so that (20) can be evaluated explicitly is another 
■attar, since the complexity of the problem is now hidden in the transformed 
fQ and V. We have examined a number of models and the most tractable requires 
that V ■ constant, in evaluating the moments u and T, (9) and (10), which are 
crucial in the aolution (20) because fQ depends on u and T. For the rest of 
this paper, we shell restrict our discussion to the "simple" model: v - 
constart* 



In this "simple" model, the explicit solution is reduced to solutions for 
u and T of the following set of coupled integral equations derived by substi- 
tution of (20) into (9) and (10): 


u(x't) e'*^ - / V f(x',v',0)/n 


+ V / dT e'^^ f d^v* V ffx(x,v,x) 
0 ^ 


( 22 ) 


[3T(x’,t)/mJ e'^‘ « / d\' Iv - 


+ V / dT e'^^ / d^v' (v - u]^ f. (x,v,t) 
0 ^ 


(23) 


As MB have stated earlier, the dependence of v and x on (x’,v',t) are sup- 
pressed in (22) and (23) for brevity. To reduce these integral equations into 
tractable form we observe that the velocity space integrals can be performed 
by applying the relations (13) - (16) to (7) end use the identity 


/♦ 

(v - u) ■ 


Ivl 


Uq sin Qt - u^ cos Qt 


4- U 2 sin Qt]‘ 


U 


+ [v^ + Ujj (1 - co« nt) - Uj^ sin flt - U2 cos + [v^ ~ (24) 


i 


where up ■ cEj^/B Is the drift speed* Carrying out the velocity space inte- 
grals of (22) and (23) is somewhat tedious; for purposes of illustration, we 
will exhibit only the procedure for (22): U3 ■ 0 and 


“1 

r*' -VT 

■ (eE /m) / dT e 

cos Qt 

.“2. 

^ 0 

-sin Qt 


( r.in at 

\ I"* 


■cos nt"| fn1 -vt 
sin Qt J l^vj ® 



(25) 


It is interesting to note that (22) is decoupled from (23) because it turns 
out that / d^v' V (x,v,t) ■ u (x',t), giving the simpler results (25). At 
t”** » uj and U2 approach the well-known Pedersen and Hall drifts respectively, 
as we %fould expect. Using (24) and (25), the integral equation (23) for the 
thermal energy 3T (x',t)/m can be %nritten explicitly as 


(3T/m) • (3 Tq/x) + 2 uJ {(1 - e"^'’^)/2 - (1 - e"^''*')/3}/(v^ + 0^) 


+ 2u^ V fi {-e sinfit + [vQ(l-e cosflt) + (2v^f Q^) e sitd^t]/ 


(4v^ + fl^))/(v^ + ft^) 


(26) 


The thermal energy of the luns is the sum of the initial thermal energy 3To/m. 
a positive monotonic contribution due to collislonal conversion of the elec- 
tric drift and a term oscillating at the ion cyclotron frequency. The oscil- 
latory term is unimportant in the F-reglon because it averages to order 


L 
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(v/n)^ which It Duch Ittt than 1* At t 

(3T/m) ♦ (3 Tq/iii) + 70 ^ Q^/(v^ + a^) (27) 

2 

Thutf tbout two thlrdt of the drift kinetic energy oiUq /2 goet into theroal 
energy* The exact amount of energy convertion depandt on model aaaumptiona; 
we regard the retult of thia model only at a guide to ion heating in the F- 
region. According to (27) then, the thermal energy increaae of the bulk of 
toptlde F-region lone under the influence of the auroral electrostatic field 
it roughly 


(T “ Tq) ^ mc^ (Ej^/B)^ ^ 8 eJ K 10“^ eV ( 28 ) 

where la in unite of mV/m. Thus, if E^^ ~ 500 mV/m [Mozer, 1981], F-reglon 
oxygen thermal energy will be ~ 2 eV - eufflclent to energize bulk oxygen 
expanaion into the topaide auroral ionoaphere [Lockwood and Titheridge, 1981], 
but not aufficient to cauae 0*^ eacape which requires > 10 eV. Note that this 
mechanism is Important only on auroral field lines where E^^ is large; if we 
consider the large-scale convection field (E^ < 100 mV/m) the thermal energy 
increaae is only < 0.1 eV « 10^ K). This result is in agreement with the 
fluid calculations of Fedder and Banks [1972]. Further, this mechanism favors 
thermal energisation of oxygen over hydrogen as observed by Ghielmettl et al. 
[1978]. 
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IV. SUPERTHERMY ANU ANISOTROPY 


Wfl presented the formulation and solution of the kinetic model in the 
previous section with a discussion of thermal energy increase of F**reglon 
ions. Since T is the thermal energy moment of the entire distribution func- 
tion, as defined in (10), the above discussion is not applicable to considera- 
tion ot non^faxwelllan features of the distribution function. These are 
discussed in this section. 

To show that the solution (20) contains a superthermal "tail" which Is 
non-Maxwellian, we need to express f(v',t) in an explicit form. This is very 
difficult, since the time convolution integral of the transformed fg is very 
complicated because u and T are complicated functions of time. Instead of 
giving a detailed numerical study, we shall show that a superthermal non- 
Maxwellian "tail" appears at early times: vt « 1 and Qt » 1. Since Ot » 
1, the oscillatory terms in u and T can be ignored because they average to 
order (v/0) or smaller. Thus, to order (v/Q), (25) and (26) yield 

Uj * c(Ej^/m)(v/Q) ; U2 « - c(E^/m); W3 ■ 0 (29) 

T/m • Tjj/m + 7 uj v^t^; u^^ - cE^/B (30) 

Defining T^/m 3 u^, we can write fQ (v',t) as 

fQ « (m/2iiTQ)^^^. (1 + Tg v^t^/TQ)"^^^ • exp - ^ (^-u)^/(l T^ v^t^/T^) (31) 

For (Tj/Tq) (vt) « 1, which is soiMwhat more restrictive than vt « 1, we can 
expand fQ above and write an approximate expression of f to order (vt)^. 
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f ■ f(t*0){« + ve j di (Tg/TQ)(vT)^Im(v' -u) - 3 Tq1) 

0 

• £(t-0){l +1 (vt)^ (Tg/Tj)I«(v*-J)^“ 3 Tq] + 0[(vt)^J} (32) 

where f(t-O) le the Mexwelllen (12). The fector in curly brackets in (32) 
indicates clearly that f ia a non'41exwellian with a auperthermal tail in the 

4 ' ^2 

superthenal regions of velocity apace where n(v'’ni) > 3T^. Mote that par** 
tides in the ai(v'**u) < 3Tq regions have been shifted to the auperthermal 
region^ 

A second important feature of the early-tiaM approximation (32) to the 
distribution function is that the supertheraal factor (curly bracket) is 
anisotropic because of (29). The anisotropic factor (v'-u)^ is thus roughly 

(^ - u)^ « [vj - c(Ej^/m)(v/Q)j2 + [V* + c(E^/n)J^ -I- (33) 

Because v/Q « 1 in the F*region, we expect the auperthermal part of the 
population in valocity space to be enhanced in the x 1 drift direction; 

A 

i.St> the 2-direction. More discussion on the anisotropy of the auperthermal 
population will be given in the next section. 

The property of anisotropy is not limited to the auperthermal part of the 
population. Indeed* wa now show that the thermal energy of the bulk of the 
distribution is also anisotropic with respect to the magnetic field. We split 
the thermal energy moment* (10) or (23)* into perpendicular and parallel 
components: 


5 / d^v' (v| - Uj^)^ f/n^ 


(34) 
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W| S / (v» - 03 )^ f/llQ 


(35) 


W i 3 T/b - + W| 


( 36 ) 


where W Is given explicitly by (26). A tedioue celculation gives 


W. - (2 - e''^)(V«) +i • U0V3 + 2 vV/(Av^-KI^)] 


[Q^ + 2 v^ (l-coeQt)] e”''^ + IQ^ - (Q cosQt + 2 v siidlt) 2v^n/(4v^-Ki^)l 


- (Q^/3) e"^''*} 


( 37 ) 


Froa W end we can obtain frcm (36). The coaplexity of (26) and (37) due 

to the rapidly oscillating teras is non-essential. We can set the sinRt and 

(1 - c>'slt) taras to xero by averaging over a tiae long coapared to cyclotron 

tlae but short coapared to collision tiae and denote t.w tiae-averaged thermal 

2 2 

aoaents with an overbar. Further, with the approxiaatlons v « Q and ITq/ih 
2 

« Uq, «m obtain' the thermal anisotropy ratio 

Wj^/(2 Wj) - 1 + I e"'’V(l - e"''*^) ; t*0 (38) 

%rhere the factor of 2 associated with coapensates for the two degrees of 
freedom for the perpendicular component. Note that (38) is valid for t / 1 /n 
froa onset because we have taken cyclotron averages. For t>0, the exact 
expressions yield Wj^/(2W|) ■ 1, as we %K>uld expect. Froa (38), we note that 
as t > «• the anisotropy again disappears. During the transition phase 
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(vt ~ 1), th« thenMl anliotropy can be quite large (~ 3 or ~ 4)» In the 
topalde F-reglon, the collision times can be tens of seconds to minutes , which 
arc of the order of the expected lifetiM oi auroral acceleration potential 
structures [Chiu and Schulz » 1978; Chiu and Cornwall, I960]. Thus, we would 
expect thermal ions in the topside auroral F>region to show preferential bulk 
heating in the perpendicular direction. 


V. DISCUSSION AIID CONCLUSIONS 


B«caat« of tho Mcottlty to mIm siapllfying ««tuaptiont in order to make 
the Bodol tractable, tie do not elaia any direct relevance of our aodei reaulta 
to obaervatlona except In a qualitative aenae* However, even at a qualitative 
level, it ia perhapa worthwhile to diacuaa the iaplicationa of our reaulta, 
which are liatad aa followat 

1. In the weakly colliaional regiaa of the topaide F-region (v « 0), a 
local ateady (100 - 1000 eV/e) can heat the Iona to a teaperature of 
aeveral electron volte, thua driving the bulk of theraal 0*^ to high altitudes 

(~ 1000 ka)* Thie Joule heating la not aufflclent to drive the bulk 0*^ to 
eacape teaperature 10 eV). 

2 

2. The therael energy gain acalea aa a %rhere a la the aaaa of the 
doainant F-reglon Ion. Thua, the aechanlaa prefera O'** heating. 

3. An Initially iaotropic Maxw<;lllan diatrlbutlon la driven to an 
laotroplc aayaptotlc atate with a higher teaperature by an external Im- 
poaed on the ayatea. In the tranaitlon atate, the dlatrlbutlon la non- 
Maxwellian with a aupertheraal tall. The Ion dlatrlbutlon In thla period (in 
the vicinity of the t « t drift energy y a u^) ia enlaotropic with reapect to 
the aagnetic field. 

4. The theraal energy increaae of the entire ion population in the 
tranaition atate ia alao aniaotropic with perpeudicular theraal energy aver- 
aging aeveral tiaea that of the parallel theraal energy. Thia aniaotropy 
diaappeara after aeveral collialon tiaea (~ ainutea), when the tranfAtlon froa 
initial to aayaptotlc etatea la coaplete. 
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Recently Ljckwood and Titherldge (1981) discussed, from the standpoint of 
observed transition altitudes, the necessity to assume an ionospheric O'*" 
heating mechanism to raise the 0*^ temperature to several eV in the auroral 
region. If Mozer's observation [1981] of large Ej^ in the topside ionosphere 
is confirmed, we believe this can be accomplished by Joule heating in the 
auroral region. We emphasize that our discussion of Joule heating (heating of 
the bulk 0*^ population) is in qualitative agreement with the calculation of 
Fedder and Banks (1972) if we apply our results to the large-scale convection 
electric field which maps to tens of mV/m in the ionosphere, rather than to 
the auroral electric field. Auroral electric fields of larger magnitude ace 
more localized and do not map through to the lower ionosphere [Rich et al., 
1981]. A probable reason for this may be the effect of the ionospheric con- 
ductivity profile on the mapping of electric fields inside the ionosphere 
[e.g. Chiu, 1974]. Joule dissipation of strong electric fields in the topside 
ionosphere, such as proposed here, may be a second reason for the absence of 
strong electric fields In the Ionosphere. 

Optical observations of the auroral 6300 A line of 0(^D) Indicate that 
the intensity Is too high (by about one order of magnitude) to be explained by 
either electron impact on atomic oxygen and/or by dissociative recombination 
of 0^ [Sharp et al., 1979]. If the bulk of the O'** ions on the topside can be 
heated to > 2 eV or if there is sufficient superthermal flux of 0*^, Increased 
population of the 0(^D) state can easily be accomplished by the charge ex- 
change interaction between the hot (> 2 eV) 0"*“ and cold atomic oxygen - the 
reaction-product atomic oxygen can easily be in the state [private communi- 
cation, A. B. Christensen]. 

Possibly our model for proc'acing superthermal 0*** by auroral electric 
fields is a first step toward a theory of auroral ion beams and conics. While 


«>• auit conatantly b« reminded of Che extreme slmpllflcetlon of the model, the 
•uperthermel part of the klnetic''rea|'onae dlatributlon function of our model 
poaaeaaea essentially the attributes of the conic distribution since the 
particlea that expand upwards along the magnetic field will acquire a 
parallel velocity from the conservation of the magnetic moment* It remains to 
be shown, however, that the model auperthermal flux at keV energies is 
sufficient to account for the observed flux. We hope to do a thorough numeri- 
cal analyaia of the diatribution in the future. Clearly our model is m^t 
intended to deal tilth high altitude (» 1000 km) phenomena wher^, we expect, 
little collisional influence and wave-particle interaction to be important. 
Observatlonally, it would be interesting to see if measurements of low-energy 
anisotropic auperthermal populations of ions at low altitudes (~ 400 km) 
(Whalen et al*, 1978] can be extended to energies as low as several eV, as In 
Klumpar (1979]. We believe, as do Lockwood and Ticheridge [1981], that n 
complete deacrlption of auroral 0^ energisation and injection into the magnet- 
osphere must begin with kinetic processes in the F-reglon ionosphere itself. 

After completion of this study, it was brought to our attention that 
similar Ideaa and formulation were diacuased in a preprint Issued by the Space 
Research Inatitute of the Soviet Academy of Sciences [Zakharov et al.,1980]. 
The discussions of the Russian work do not involve solutions of the Krook 
model as is done here* 
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APPENDIX A: TRANSFORMATION OF THE BOLTZMANN OPERATOR 

Application of the transformation (x,v,t) (x',v',t), (13) - (16), to 

A 

the components of the Boltzmann operator where the curly brackets 

indicate functionals, yields the sought -for identity 




(A-U 


provided the force is given by 





+ e(v X Bz)/c 


(A-2) 


as assumed in our model* 

The calculation is straightforward hut somewhat tedious, so we shall only 
provide the salient points here* Direct differentiation gives 


3f (x,v,t) 
■St 


V £(;^«,^’,t) + . V^, f(^,v',t)(A-3) 


where 7* and 7^, are gradient operators with respect to x* and v' respective- 
ly* From (13) - (16), we have 


3x' 

dt 


v(x',v',t) 


(A-4) 


3v’ 

IT 


l(^) 


(A-5) 


Next, 
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V • V£(x.v,t) - 


(A-6) 


1.1 

b«cxuse • Application of (13) - (16) once aoro gives 

^(v) • f(x,v,t) - f(v’) • V^, f(x',v',t) (A-7) 

The identity (A-1) is obtained by using (A-3) - (A-7), 
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APPENDIX B: CONSERVATION OF PARTICLES 


It ii pointed out In Section III that the choice of model initial and 
asymptotic distribution functions^ (12) and (7) respectively! nuat be self- 
consistent with the imposed model constraint V • - 0 which amounts to (21) 
via (4): 


n(x',t) - / d\' f(J',^,t) - Uq (B-l) 

where f is given in terms of the forms (7) and (12) by (20). Here we prove 
(B-l) by direct integration. 

Substitution of (20) into (B**!) yields two Integrations over d^v*. The 
velocity space integration over f(x'!V*!0)! (12), yields straightforwardly 

nQ. The velocity space integration over Iq looks very complicated but inspec- 
tion of (24) shows that 


(V - $)^ - l^ - w(t)l^ (B-2) 

where w(t) does not depend on v*. The Maxwellian form of fg now allows 
/ d v' fp to be integrated! yielding a result Independent of w. Thus, (B-1) 
yields 

_ t 

n(x*!t) - iIq e"'’^ + "O '' " Sq 
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